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ABSTRACT 

We report on evidence for orbital phase-dependence of the 7-ray emission 
from PSR B 1957+20 black widow system by using the data of the Fermi Large 
Area Telescope. We divide an orbital cycle into two regions: a region containing 
the inferior conjunction, and the other region containing rest of the orbital cycle. 
We show that the observed spectra for the different orbital regions are fitted 
by different functional forms. The spectrum of the orbital region containing 
inferior conjunction can be described by a power-law with an exponential cutoff 
(PLE) model, which gives the best-fit model for the orbital phase that does not 
contain the inferior conjunction, plus an extra component above ~ 2.7 GeV. The 
emission above 3 GeV in this region is detected with a ~ 7a confidence level. 
The 7-ray data above ~ 2.7 GeV are observed to be modulated at the orbital 
period at the ~ 2. 3d level. We anticipate that the PLE component dominating 
below ~ 2.7 GeV originates from the pulsar magnetosphere. We also show that 
the inverse-Compton scattering of the thermal radiation of the companion star 
off a "cold" ultra-relativistic pulsar wind can explain the extra component above 
~2.7 GeV. The black widow pulsar PSR B 1957+20 may be the member of a new 
class of object, in the sense that the system is showing 7-ray emission with both 
magnetospheric and pulsar wind origins. 

Subject headings: binaries: close; — magnetic fields: pulsars; — stars: neutron; 
- pulsars: individual (PSR B1957+20) 
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Introduction 



The Fermi Large Area Telescope (Fermi) has me asured the pulsed 7 -r ay emission from 
millise cond pulsars, which is magnetospheric in origin ( Abdo et al. 2009al lb1; IGuillemot et al. 
2012aJ). Furthermore, more millisecon d pulsars (MSPs) have been discovered as Fermi 



unidentified sources in radio bands (e.g.JCognard et al.ll201ll ; iRansom et al .1120 lit IGuillemot et al 
2012a ; Kerr et al. 2012 ; Ray et al. 2012b . suggesting MSPs are one of the major Galactic 
sources in the high-energy 7-ray sky. The pulsed 7- ray emissions from energetic MSP, PSR 
B 1957+20, was reported bv IGuillemot et all (boi2b| ). PSR B 1957+20 has been known as the 
first "black wido w" binary system, in which the MSP is destroying the low mass companion 
star (< O.1M ) (IFruchter. Stinebring &; Taylorlll988l ).So far, at least 10 black widow syst ems 
have been discovered in the Fermi unidentified sources jRobertsll20l"ll ; iRav et alllioii l. 



Stappers et a 



(|2003f) detected the unresolved X-ray emission around PSR B1957+20, 



and lHuang et al.l (120121 ) found the modulation of the unresolved X-ray e missions with the or- 
bital 



bital period of P n ~ 9.2 hrs, suggesting the origin of an int ra-binary shock (iHarding fc Gaisser 
1990 ; lArons fc Tavani 1993 ; Cheng. Taam fc Wang![2006 ). In the intra-binary shock scenario 



the matter of the wind from the extremely low mass companion is injected by the irradiation 
of the strong wind /radiation from the M SP, and collides with the pulsar wind at a position 
not far from the companion (IChengill989l ). The pulsed radio emissions from PSR B1957+20 
show s the eclipse lasting for about 10% of the orbital phase (IFruchter. Stinebring &; Taylor 
19881 ). The corresponding size of the eclipse is on order of ~ 10 11 cm, which is larger than 
the Roche lobe size ~ 2 x 10 10 cm , which implies that s tellar material is expelled from 
the companion by the pulsar wind ( IRuderman et al.l Il989l ) . The irradiation of the pulsar 
wind/radiation on the companion star of PSR B1957+20 system has also been expected 
by t he observations of the orbital modulation of the optical emissio ns from the compan- 
ion (jCallanan. van Paradijs fc Rengelinkl Il995l ; iReynolds et al.l 120071 1 . In addition to the 
PSR B 1957+20 system, the orbital modulation of the optical emission associated with 
the irradiation of the M SP has been confirmed for seve ral binary systems, for example, 
the J1023+0038 system ( Thorstensen fc Armstrong! 2005h. the bla ck widow system candi- 



( IKataoka et al 



2012 



Romani 



dates 2FGL J23 39.6-0 532 (iRomani fc Shawll201ll ; lKong et alJl2012ah and 2FGL J1311. 7-3429 



J1808.4-3658 faurderi et al.lhooi iDelove et al.ll2008h . 



2012|). and the accreti ng MSPs in quiescent state, e.g. SAX 



The high-energy emission associated with an intra-binary shock between a pulsar and 
a companion star has been detected from the so-called 7-ray binary PSR B1259-63/LS 2883 
system, which is composed of a canonical pulsar with a period of P ~ 48 ms and a high-mass 
Be star ( iJohnston et al.lll992l ; lAharonian et al.ll2009t lUchiyama et al.ll2009l ). In the 7-ray bi- 
nary, it has been proposed that the shock stands at the interface between the pulsar wind and 
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the Be wind/disk, and un-pulsed radio to TeV radiations are produced via the synchrotron 



1997 




Dubus 


2006; 


Takata et al. 


2012) 



Takata fc Taaml 120091 ; Kong et all 1201 ll ; Kong. Cheng fc Huand 12012b 



vation during the 2010/2011 periastron passage (lAbdo et al.l 1201 It iTam et al.l l201ll ). Dif- 



ferent models have been proposed for explaining the flar e-like GeV emission from the P SR 



B1259-63/LS 2883 system: the Doppler boosting model (IKong. Cheng fc Huang! l2012bl ) as 



sumes that the flow in the bow-shock tail is r elativistic and the synch rotron emission is 



Doppler-boosted, while the pulsar wind model (IKhangulyan et al.l 120121 ) assumes that the 
GeV emissions are produced by the inverse- Compton process of the cold-relativistic pulsar 
wind. 

The 7-ray emissions associated with the pulsar wind in the black widow system PSR 
B1957+20 have not been reported yet. However, we note that for a typical 7-ray binary 
system, the distance to the shock from the pulsar is order of R s ~ 0.1 — 1 AU, corresponding 
to R s /R\ c ~ 10 4 - 10 5 , where R lc = 2.4 x 10 8 (P/0.05 s) cm is the light cylinder radius 
and P is the spin period of the pulsar. For the black widow system, the shock stands at 
the distance R s ~ 10 10 — 10 11 cm, corresponding to R s /R\ c ~ 10 3 — 10 4 , which is roughly 
one order of magnitude smaller than that of the 7-ray binaries. Therefore, the black widow 
system will provide us a unique laboratory to probe the physics of the pulsar wind, and may 
be candidate of 7-ray emissions from the pulsar wind. In this paper, we report searching for 
orbital modulation of GeV 7-rays from PSR B1957+20 system in the Fermi data to obtain 
an evidence of the 7-ray emissions from the intra-binary space. 



2. Data analysis and results 

In this analysis, Fermi-LAT data were taken between 2008 August 4 and 2011 De- 
cember 5. We restricted the events in the "Source" class (i.e., event class 2) under the 
"P7SOURCE_V6" instrumental response functions (IRFs). In addition, we excluded the 
events with zenith angles greater than 100° to reduce the contamination by Earth albedo 
gamma-rays. The data were analyzed using the Fermi Science Tools v9r23pl, available 
from the Fermi Science Support Center 0. Events were selected within a circular region 
of interest (ROD centered at the posi tion of: R.A. = 19 h 59 m 36!77, decl. = 20°48'15'. / 1 



(lArzoumanian. Fruchter fc Taylorl Il994l ). In order to reduce systematic uncertainties due 



to the surrounding complex galactic region and to achieve better background modelling, 



1 http : / /fermi . gsfc . nasa.gov/ssc /data/ analysis / software / 
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photon energies were restricted to be above 200 MeV and a size of radius 5° of the ROI 
was adopted throughout the analysis. As an attempt to search for evidence of orbital mod- 
ulation, an orbital phase wa s assign ed to each 7-ray photon based on the timing ephemeris 
reported by Guillemot et al. ( 2012b ). using the Fermi plug-in for Tempo2§[ taking an aper- 
ture radius of 1°. No sign ificant evidence of orbital m odulation at a > 2a level was found 
by e mploying the H-test (Ide Jager &: Buschingl 120101 ). which is consistent with the results 
from iGuillemot et al.l (|2012bJ). To access the effect of the sky exposure as a function of time 
throughout the observation, a lightcurve was created with a bin size equal to the one-tenth of 
the orbital period of the binary system, the exposure was calculated with the tool gtexposure. 
The summed exposure in each bin was found to deviate from the mean value by less than 
2%, suggesting that the variation in exposure has a minute effect on the orbital lightcurve. 
To avoid bias due to the alignment of the first bin with the orbital phase zero, the bins were 
shifted by phases of 0.2, 0.4, 0.6 and 0.8. The results after shifting the lightcurve do not 
change our conclusion. 



Recently, iHuang et al.l ( 120121 ) have discovered orbital-modulated X-ray emission from 
the PSR B1957+20 system. To investigate whether the 7-ray spectral properties also vary 
with the orbital phase ((f)), we divided an orbital cycle into two parts: half of the orbit 
centered at the superior conjunction (hereafter "Phase 1") and the other half containing the 
inferior conjunction (hereafter "Phase 2"). Unbinned likelihood analysis was performed for 
both regions with the help of gtlike. A spectr al-spatial model co ntaining all other sources 
reported in the Second Fermi source catalog ( INolan et al.l 12012k 2FGL hereafter), within 
10° from the center of the ROI was used for source modelling, resulting in a total of 13 
point sources in the model. All the 2FGL sources were assumed to have the spectral type 
as suggested in the catalog. The spectral parameters of all sources within the ROI were 
set to be free, while the parameters for sources outside the region were kept fixed. We 
also included the Galactic diffuse model (gal_2yearp7v6_v0 . f its) as well as the isotropic 
spectral template (iso_p7v6source . txt) and allow their normalizations to be free. Since the 
7-ray spectra of pulsars in general can be characterized by a powe r-law with an exp onential 



cut-off (PLE) model, dN/dE = N (E /0.2 GeV)" r exp [-(E/E c )] ( lAbdo et al.ll2010j ). where 



N is a normalization factor, T is the photon index and E c is the cutoff energy, we performed 
fitting of the spectrum average over each orbit segment with a PLE model. For Phase 1, 
the best-fit parameters are a photon index of Ti = 0.83 ± 0.56 and a cut-off energy of 
£ C utoff,i = 0.84 ± 0.30 GeV, while for Phase 2, T 2 = 1.58 ± 0.30 and £ cut off,2 = 2.28 ± 0.88 
GeV. A summary of the fitting results is presented in Table [TJ The difference in cut-off 
energies in the two phases suggests that the spectrum of Phase 2 extends to higher energies 



See http: / / fermi . gsfc . nasa. gov /ssc/ data/ analysis / user /Fermi_plug_doc .pdf 
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than that of Phase 1. Moreover, the observed spectrum of Phase 1 is fitted by the PLE model 
with a Test-Statistic (TS) value of 205, while it is fitted by a PL model with a TS value of 
168. This indicates that the PLE model is more favored over the PL model in describing 
the spectrum of Phase 1 at a ~ 6.1a confidence level. For Phase 2, the TS values of the 
best-fit PLE and PL models are 161 and 143 respectively, indicating that the PLE model 
is not significantly favored over the PL model. The spectral energy distributions (SEDs) 
for both orbital phases are presented in Figure [TJ We computed the extrapolated fluxes in 
the energy range 0.1 - 300 GeV for both Phase 1 and Phase 2, which ar e also presented 



i n Tabl e [TJ The values of the fluxes are consistent with that reported in iGuillemot et al. 



(l2012bl ) within uncertainties. We note that in the source model, the best-fit log-parabola 
model for the source 2FGL J1949. 7+2405 yields a curvature coefficient (3 « 0.2 - 0.3, which 
is different from th at given in th e 2FG L catalog, where the value of j3 is fixed at 1. However, 



it is pointed out in iNolan et al.l ( 120121 ) that such highly curved spectrum is not necessarily 
robust, possibly due to the densely populated region near the Galactic plane. Moreover, a 
likelihood analysis using an ROI of radius 15 deg (see below) yields a similar result, suggesting 
that the deviation could be justifiable. 

Since the radius of the point-spread function (PSF) is relatively large near the low 
energy bound of 200 MeV, the spectral results may be affected by the wings of the PSFs of 
surrounding sources. Thus, we performed binned likelihood analysis with an ROI of radius 
15°around PSR B 1957+20. The spectral fitting results are presented in Table [TJ These 
results are consistent with those obtained using an ROI of radius 5°. We also note that the 
difference in the best-fit cut-off energies in Phase 1 and Phase 2 agrees with the conclusion 
on the difference in the spectral properties of the two orbital phases. 

By comparing the 7-ray spectra, we speculate that the spectrum for Phase 2 could be de- 
scr ibed by two components , one being magnetospheric, which has also been well-established 



by IGuillemot et al.l (j2012bl ). and the other component above ~ 3 GeV, coming from the 
interaction between the pulsar wind and the companion star. To describe this component 
while avoiding bias towards any emission scenario, we adopted a simple Gaussian profile 
dN/dE = Aexp — [(E — E) 2 /<Jq] to fit the data along with a PLE model with the spectral 
index and the cutoff energy fixed at the best-fit values derived from Phase 1. The best-fit 
two-component model is presented as the solid line in the right panel of Figure [TJ The best-fit 
parameters for the Gaussian component are E = 3.76 + 0.59 GeV and oq = 1.10 + 0.39 GeV. 
To access the significance of emission of this extra component, we performed likelihood anal- 
ysis using the data at energies above E — <jg ~ 2.7 GeV by assuming a simple PL model, 
alleviating the effect of the spectral model on the significance. For Phase 2, a TS value of 
55 was achieved, corresponding to a significance of ~ 7a, whereas for Phase 1, a TS value 
of 14 was obtained, indicating that the detection significance is below 4a at > 2.7 GeV for 
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Phase 1. We computed the TS maps using gttsmap at energies > 2.7 GeV for both Phase 1 
and Phase 2, which are shown in Figure [2](a). The difference in the TS values for the two 
orbital phases is readily observable by comparing the maps. 

In addition, we investigate the relation between the rotation phase of the pulsar and 
the variation of spectral characteristics in the two orbital phases by removing the pulsed 
component. A pulse profile of PSR B1957+20 was constructed by phase-folding the photons 
with energies above 0.2 GeV and within an angular distance of 1° from the position of the 
pulsar. The gamma-ray peaks in the profile was then fitted with two Lorentzian functions. 
Data were excluded within the full widths at half-maximum of the fitted peaks, which consist 
of phases 0.103 - 0.215 and 0.598 - 0.620. The pulse profile and the fitted Lorentzian 
functions are illustrated in Figure El We compare the significance of emission in Phase 2 
below and above 2.7 GeV, before and after removing the pulsation peaks. At energies below 
2.7 GeV, the TS value decreased from 105 to 38, while at energies above 2.7 GeV, the TS 
value decreased from 55 to 36. Figure EJ^b) and (c) show the TS maps for a visual comparison 
before and after the removal of the pulsed component. 

The above results indicate the presence of emission above 2.7 GeV, which is significantly 
detected in Phase 2 but not in Phase 1. Hence, the 7-ray emission from PSR B 1957+20 is 
dependent on the orbital phase of the binary system. Moreover, the removal of the pulsed 
emission component has caused a greater decrease in the detection significance at energies 
below 2.7 GeV than above, implying the majority of the emission of 2.7 GeV in Phase 2 are 
not produced inside the pulsar magnetosphere, but in the intra-binary region. 

Since the best-fit two-component model suggests that photons at > 2.7 GeV may be 
modulated with the orbital phase, we constructed a phase-folded lightcurve and derived the 
significance of the modulation. To allow for enough photon statistics we selected events 
with energies greater than 2.7 GeV, the size r of the aperture. The best profile (r = 0?965 ) 



is presented in Figure E2 with an H-test Test Statistic of 19 ( Ide Jager fc Buschingl |2010| ). 
corresponding to a significance of ~ 3.3a. The post-trial significance, associated with 21 
trials on the aperture radius, is ~ 2.3o\ 



3. Discussion 

We have reported the evidence of orbital-modulated GeV 7-ray emissions from the black 
widow system PSR B1957+20. We have fitted the spectrum of Phase 2 (half of the orbit 
centered at the inferior conjunction) using best-fit PLE model of Phase 1 plus the Gaussian 
component with E = 3.76±0.59 GeV and ac = 1.10±0.39 GeV. The additional component 
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in the spectrum has been found with ~ 7a confidence level. We have also shown that the 
orbital modulation above ~ 3 GeV can be found with a ~ 2.3a confidence level. 

It can be expected that the PLE component below 3 GeV will be dominated by the 
emissions from the pulsar magnetosphere. For the additional component above 3 GeV in 
Phase 2, we will speculate the possibility of the inverse-Compton process of a "cold" ultra- 
relativistic pulsar wind. In this section, we will calculate the expected spectra of the emissions 
from pulsar magnetosphere and from the pulsar wind. 



3.1. Outer gap emissions 



For the magnetospheric emissio n, we apply the two-layer outer gap model explored by 
Wang. Takata fe Cheng! (120101 . 1201 ll ) . The dashed lines in Figure [5] represent the calculated 
phase-averaged spectra of the curvature radiation from the outer gap. The results are for 
the inclination angle of a = 60° and the viewing angle of £ = 80°, respectively, which 
were assumed to explain the pulse profile and the spectrum. A larger viewing angle is 
required to explain the phase separation (~ 0.45 orbital phase) of the two peaks in the 
light cu rve of the pulsed emissio ns. The assumed viewing geometry is consistent with that 
used in I Guillemot et al.l ( 12012b! ) . For the present outer gap model, the 7-ray luminosity 
is expected to be order of L gap ~ f 3 L sd , where / is fractional gap thickness, which is 
ratio of the gap th ickness to the light cyl inder radius at the light cylinder, and L s< i 



( Manchester et al. 2005 ) 3 ] is the spin down power of PSR B1957+20. 



1.6 x 10 35 erg s 1 

We no te that because t he period time derivatives before and after correction of Shklovskii 

-20 



effect (phklovskii 



( Manchester et al. 



9701) due to proper motion are P ~ 1.68 x 10 20 and ~ 1.13 x 10 



20051 ). respectively, the spin down power and hence our main conclusion 



are not changed much by the Shklovskii effect. For the fractional gap thickness, we apply 
the outer gap model controlled by the magnetic pair-creation process near the stellar surface 
dTakata. Wang fc CheneJboiof ). which implies / ~ 2.5(P/0.1s) 1 / 2 ~ 0.3 for PSR B1957+20. 
The right panel in Figure compares the model spectrum of the outer gap and the observed 
7-ray spectrum in off-peak orbital phase in Figure HJ where we expect the magnetospheric 
emission dominates. We can see in Figure |5] that the emissions from the outer gap can 
explain well the Fermi data in off-peak orbital phase. 



3 http://www. atnf.csiro.au/research/pulsar/psrcat 
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3.2. Pulsar wind emissions 



It has been thought that a "cold" pulsar wind, which will be mainly composed of the elec- 
trons and positron pairs, is injected into the interstellar space by the pulsa rs. The pulsar wind 
is acc elerated beyond the light cylinder to a Lorentz factor Fyy ~ 10 5-7 (IKennel fc Coroniti 
19841 ) . If the cold ultra-relativistic pulsar wind interacts with the inter-stellar medium, the ki- 
netic energy of the pulsar wind is converted into the internal energy at the shock. The shocked 
pulsar wind can emit X-rays via synchrotron radiation and is observed as the pulsar wind 
nebula. The un-shocked cold pulsar wind is dark in the X-ray bands, because the wind does 
not emit synchrotron photons. However, it can been thought that the cold ultra-relativistic 
pulsar wind emits very high-energy 7-rays by the inverse-Compton process between the cold 
pulsar wind and the background soft-photons from the pulsar magnetosphere or from the 
companion in the binary system. In fact, studies hav e been carried out on the 7-ray emis- 
sions from the cold pulsar wind of the Crab pulsar (lAharonian. Bogovalov fc Khangulvan 



20121) and of the 7-ray b inary PSR B1259-63 ( IKhangulyan et al 
Khangulvan et al.ll2012[ ). 



2007 



Takata fc Taam 



2009 



For PSR B1957+20 black wido w system, the temper ature of the heating side of com- 
panion is as high as T c s ~ 8300 K ([Reynolds et al.l 120071 ) . Roughly speaking the luminos- 
ity of the inverse-Compton process may be estimated by L IC ~ (1 — cos 9 )aTn p ha s L sd ~ 
5 x 10 _3 (1 — cos6* )(T cff 3 /8.3) 3 i?^ 10 a~]; 1 L sd , where 9 is the collision angle between the pulsar 
wind and the soft photon, R s ~ 10 10 i? s> io cm is the Roche lobe radius and a s = 10 11 a s .n cm 
is the separation between the pulsar and the companion. In addition, n p h ~ a s bT^ s B? s / kbca 2 
with Izb being the Boltzmann constant, is the typical number density of the soft photon field, 
and T e ff 3 = (T ef f/1000 K). Because the typical luminosity of outer gap emissions is of the 
order of L gap ~ f 3 L sd ~ 0.03L sd , the inverse-Compton process of pulsar wind may produce 
7-ray radiation with a flux level that is of the same order of magnitude as that of the gap 
emissions, if the inverse-Compton process occurs with head-on collision, which can happen 
around the inferior conjunction, because the companion is located between the pulsar and 
the Earth. In fact, we have measured the significant excess of the 7-ray emissions above 3 
GeV during the inferior conjunction passage. The typical energy of the scattered photons is 

eff 



3r^fcT eff ~ 5 GeV(T eff , 3 /8.3)(iV/5 x 10 4 ) 2 . 



We perform a more detail calculation for the inverse-Compton process between the cold 
ultra-relativistic pulsar wind and the soft-photon field from the c ompanion. The emissivity 
of the inverse-Com pton process of a particle may be expressed by (IBegelman fc Sikoralll987l ; 
Takata et al~ll2012h . 



dPic 



V 2 



f3 cos 9 )Ib/h 



KN 



dfl , 
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where da' KN /dQ' is the differential Klein-Nishina cross section, f3 = yT^T--T/IV, V = 
r^(l — /3cos#i) _1 with 9i being the angle between the direction of the particle motion 
and the propagating direction of the scattered photon, h is the Planck constant, h is the 
background photon field and 8 C expresses the angular size of star as seen from the emission 
point. 

We assume that all pulsar wind's particles have a Lorentz factor of IV = 4 x 10 4 . We 
consider that the energy flux of the p ulsar wind depends on the colatitude (9) measured from 



the rotation axis as L w (6) oc sin 9 (IBogovalov k, Khangoulianll2002l ; iLyubarskyl l2002f ) . and 



that total power of the wind corresponds to the spin down power. We also assume that half of 
the energy flux is carried by the kinetic energy of the pulsar wind. The fitting of the optical 
modulation using ELC model implies that the inc lination angle of the orbital plane with 



respect to the sky is i ~ 65° ([Reynolds et al.l 120071 ). In the present calculation, we assume 



i = 67°, and the viewing angle £ = 80° respective to the rotati on axis of PSR B 1957+20. The 



distance to the source has be en estimated as<i~2.5±lkpc (IFruchter. Stinebring fc Taylor 



19881 ; iGuillemot et al.ll2012bl ). In the present calculation, we apply d = 2 kpc. We note that 
the observed flux can be reproduced more easily with a larger inclination angle (i) and a 
smaller distance than the estimated values within the tolerance of their uncertainties. 

The spectra of the inverse- Compton process of the pulsar wind are described by the dot- 
ted lines in Figure [5j and the temporal behavior is summarized in Figure El where the phase 
zero corresponds to the inferior conjunction. We have expected that the pulsar passes the 
inferior conjunction at the center of X-ray peak phase (Phase 2). At the inferior conjunction 
passage, because the companion is located between the pulsar and the Earth, the inverse- 
Compton process occurs with head-on collision rather than tail-on collision. We found that 
the emissions from inverse-Compton process dominates the curvature radiation from the 
outer gap in the spectrum above 3 GeV, and can explain the observed excess in the Fermi 
data at the peak orbital phase, as the left panel in Figure \5\ indicates. In the off-peak phase 
(Phase 1), because the inverse-Compton process occurs with tail-on collision rather than 
head-on, the inverse-Compton process is suppressed. Hence, only magnetospheric emissions 
contribute to the observed emissions, as shown in right panel of Figure [5j 

As we have seen, the inverse-Compton process of the cold ultra-relativistic pulsar wind 
can explain the observed extra component above 3 GeV in the orbital phases around the 
inferior conjunction. Hence, we would like to emphasize that the PSR B1957+20 binary 
system is (1) a plausible candidate manifesting the emissions from the cold ultra-relativistic 
wind, and (2) a candidate of the new class of 7-ray binary, in the sense that the 7-ray binary 
is composed of a MSP and a low mass star, and showing both magnetospheric and pulsar 
wind emissions. 
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It is unlikely that the extra component above 3 GeV during the inferior conjunction pas- 
sage is originated from the inverse- Compton process of shocked pulsar wind. We may estimate 
the magnetic field of t he pulsar wind after the shock by B = 3L^ 2 a s c _1 / 2 cr^ 2 (l + aw)^ 1 ^ 2 



( iTakata fc Taamll2009l ). where aw is the ratio of the magnetic energy to the particle energy at 
the shock. The energy density of the thermal radiation from the companion is at most U p h = 
<7sB?cfr/ c - Hence, the ratio of powers of the synchrotron radiation and the inverse-Compton 
process of the shocked pulsar wind is estimated as Us/Uph ~ 10Lgj 2 35 (T cffi3 /8.3)~ 4 criy(l + 
aw), where L s d,35 = (L s d,35/10 35 erg s" 1 ). On the other hand, the ratio of the observed 
fluxes of the orbital modulated X-ray emissions (~ 10~ 13 erg cm -2 s -1 , Huang et al. 2012) 
and GeV emissions (~ 10 -11 erg cm -2 s _1 ) are ~ 10 -2 . Hence, the inverse-Compton process 
of the shocked pulsar wind will not explain the emissions above ~ 3 GeV observed at the 
inferior conjunction passage. 

Finally, we remark that the black widow pulsars discovered in Fermi unidentified sources 
are good targets for searching orbital phase-dependent 7-ray spectra. In particular, the black 
widow candidates 2FGL J2339.6-0532 and 2FGL J1311. 7-3429 show orbital modulation in 
optical, and it is evident that the former candidate shows X-ray modulation as well. Hence, 
they are an excellent probe for 7-ray emission with both magnetospheric and pulsar wind 
components. 

In summary, despite the limitation from the significance of each individual approach, 
the pieces of evidence as a whole has allowed us to conclude the detection of orbital modu- 
lated 7-rays from the PSR B1957+20 system in Fermi-LAT data. We have shown that the 
significant emissions above 3 GeV appears around inferior conjunction, while the emissions 
below 3 GeV are steady and are described by the pulsar emissions. We expect that the 
modulated emissions above 3 GeV is originated from the inverse-Compton scattering of the 
thermal radiation of the companion star off the cold ultra-relativistic pulsar wind. 
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Table 1. Summary of spectral fitting results of the two orbital phases 



Orbital 


ROI radius (°) 


Spectral 


r 


E c 


E 




Flux a 


phase 


(analysis type) 


model 




(GeV) 


(GeV) 


(GeV) 


(KT 9 ph cm" 2 s" 1 ) 


0.5 - 1.0 


5 (unbinned) 


PLE 


0.83 ±0.56 


0.84 ±0.30 






8.95 ± 1.79 


(Phase 1) 














(12.8 ±4.42) 




15 (binned) 


PLE 


1.21 ±0.57 


1.04 ±0.46 






8.76 ± 1.99 
(14.3 ±5.94) 


0.0 - 0.5 


5 (unbinned) 


PLE 


1.58 ±0.30 


2.28 ± 0.88 






8.47 ± 1.57 


(Phase 2) 














(15.0 ±4.52) 




5 (unbinned) 


PLE + Gaussian 


0.83 b 


0.84 b 


3.76 ±0.59 


1.10 ±0.39 


7.19 ±0.85 
(10.2 ± 1.22) 




15 (binned) 


PLE 


1.88 ±0.19 


3.17 ± 1.15 






8.12 ± 1.28 
(16.6 ± 3.77) 




15 (binned) 


PLE + Gaussian 


1.21 b 


1.04 b 


3.80 ± 0.70 


1.00 ±0.54 


6.76 ± 0.96 



(10.9±1.57) 



a Integrated photon flux in the energy range 0.2 - 300 GeV. The values in brackets represent the 0.1 - 300 GeV 
extrapolated photon flux. 

b Model parameters without quoted errors are fixed at the values given. 




Energy (MeV) Energy (MeV) 



Fig. 1. — Spectral energy distributions of 7-ray emission from PSR B1957+20. Data points 
were derived from likelihood fitting of individual energy bins, in which a simple PL is used to 
model the data. 90% upper limits were calculated for any energy bin in which the detection 
significance is lower than 3cr. Left: Spectrum averaged over Phase 1. The solid line shows 
the best-fit PLE model from fitting the data above 0.2 GeV. Right: Spectrum averaged over 
Phase 2. The solid line represents the fitted two-component model, with the PLE component 
shown as a dashed line and the Gaussian component shown as a dash-dotted line. 
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Fig. 2. — Test-statistic (TS) maps of 2° x2° regions centered at the position of PSR B 1957+20 
(labeled by green crosses). The color scale below each pair of images is used to indicate the 
TS values, (a): (i) TS map at energy > 2.7 GeV using only photons in Phase 1. (ii) Same 
as (a)(i) but using only photons in Phase 2. (b): (i) TS map at < 2.7 GeV for Phase 2. 
(ii) Same as (b)(i) but with data within full width at half maximum of the pulsation peaks 
removed (see text), (c): Same as (b) but with energy > 2.7 GeV. 
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Fig. 3. — Pulse profile of the 7-ray emission from PSR B1957+20, produced by using data 
above 0.2 GeV and within 1° from the position of the pulsar. The red solid line represents 
the fitted Lorentzian functions. The shaded regions indicate the pulse phases within the full 
widths at half maximum of the two peaks. 
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Fig. 4. — 7-ray lightcurve of PSR B1957+20 folded at the orbital period using the Fermi 
plug-in for Tempo2, with the optimized size of aperture of 0?965. Two orbits are shown for 
clarity. The shaded regions correspond to the phase of radio eclipse, which is the center of 
Phase 2. 
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Fig. 5. — The spectrum averaged over Phase 1 (left panel) and Phase 2 (right panel). The 
dashed and the dotted lines are the spectra of the emissions from the outer gap and the cold 
ultra-relativistic pulsar wind, respectively. The solid line represents the total contribution 
from the components. 
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Fig. 6. — Calculated orbital modulation for 7-ray emissions from the inverse- Compton pro- 
cess of a cold ultra-relativistic pulsar wind. The phase corresponds to the inferior conjunc- 
tion. 



